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ABSTRACT 

The comoving mass density of massive black hole (MBH) remnants from pre-galactic star formation 
could have been similar in magnitude to the mass-density of supermassive black holes (SMBHs) in the 
present-day universe. We show that the fraction of MBHs that coalesce during the assembly of SMBHs 
can be extracted from the rate of ring-down gravitational waves that are detectable by Advanced LIGO. 
Based on the SMBH formation history inferred from the evolution of the quasar luminosity function, we 
show that an observed event rate of 1 yr _1 will constrain the SMBH mass fraction that was contributed 
by MBHs coalescence down to a level of ~ 1CT 6 for 20A1© MBH remnants (or ~ 1(T 4 for 26OM 
remnants). 

Subject headings: black hole physics - cosmology: theory - galaxies: formation 



1. INTRODUCTION 

The first episode of star formation (Pop-Ill) occurred 
in metal-free gas that cooled to low temperatures through 
collisional excitation of molecular hydrogen, H2 (Bromm 
& Larson 2003). The first stars could have formed inside 
mini-halos of masses M > 10 5 M Q at redshifts z > 20 
(Barkana & Loeb 2001). Numerical simulations indicate 
that the first stars were probably more massive than the 
sun by two orders of magnitude (Bromm, Coppi & Lar- 
son 1999; Abel, Bryan & Norman 2000). Pop-Ill stars with 
a mass M* > 26OM or 25 < (M*/M ) < 140 end their 
life by forming a MBH remnant (Heger & Woosley 2002; 
Heger et al. 2003). 

The formation of Pop-Ill stars must have been self- 
regulating, as the UV radiation emitted by these stars 
dissociated H2 throughout the universe (Haiman, Rees & 
Loeb 1997) and boiled (via photo-ionization heating) gas 
out of shallow potential wells with a virial temperature 
< 10 4 K (Barkana & Loeb 1999), whereas supernovae 
explosions expelled gas from even deeper potential wells 
(Yoshida et al. 2003). Moreover, Pop- III supernovae en- 
riched star-forming regions with metals above the critical 
value of ~ 0.1% of the solar metallicity, allowing fragmen- 
tation of the gas into lower-mass stars (Bromm & Loeb 
2003). Even with this self-regulation included, Madau & 
Rees (2001) estimate that the co-moving mass density in 
MBHs at z ~ 20 may have been comparable to the ob- 
served mass density of supermassive black holes (SMBHs) 
in the present-day universe. Thus, in principle, the coa- 
lescence of MBH remnants could have made a significant 
contribution to the SMBH mass budget. 

More recent semi-analytic calculations of the merger his- 
tory of SMBHs (Islam, Taylor & Silk 2002) find that al- 
though the mass density in Pop-Ill MBH remnants may 
be sufficient to provide the present day density of SMBHs, 
most of the MBHs end up in satellites and populate galaxy 
halos. Nevertheless, hierarchical merging inevitably leads 
to the formation of early dwarf galaxies, which could sink 



via dynamical friction to make dense clusters of MBHs at 
the centers of massive galaxies that are assembled later 
(Madau & Rees 2001). These central MBH clusters pro- 
vide sites for MBH coalescence and could produce seeds for 
SMBH growth (analogous processes take place within the 
segregated cores of globular star clusters; see Portegies- 
Zwart & McMillan 2000). 

The rise of the SMBH population is traced by the 
evolution of the quasar luminosity function (e.g. Yu & 
Tremaine 2002). While the gas accretion history of the 
quasar population accounts for most of the locally ob- 
served SMBH mass, coalesced MBHs may also contribute 
some small fraction, /mbh, of that mass. For a given value 
of /mbh, the coalescence rate of MBHs corresponding to 
the buildup of mass in SMBHs may then be calculated 
based on the observed quasar history. In this Letter we 
show that detection of gravitational waves from MBH coa- 
lescences by Advanced LIGCP (LIGO-II) will constrain the 
fraction of the SMBH mass that was assembled through 
coalescence of MBHs and provide a direct probe of the 
early population of very massive stars. 

2. RING-DOWN GRAVITATIONAL WAVES FROM MBHS AND 
DETECTION BY LIGO-II 

First, we consider the expected gravitational radiation 
signal from merging MBHs at high redshifts. The coa- 
lescence process of two black-holes can be separated into 
three phases: an inspiral phase, a coalescence phase, and 
finally a ringing phase of the merged black-hole product 
(Flanagan & Hughes 1998) as it settles to a Kerr metric. 
Of these, the inspiral and ringing phases generate well un- 
derstood wave-forms. Inspiraling MBHs at high redshift 
will be most luminous at frequencies near 1 Hz, just be- 
low the sensitivity band for LIGO-II. However, following 
a merger the MBHs would ring at higher frequencies that 
are within the LIGO-II sensitivity band. 

To compute the gravitational wave strain for black-hole 
ringing we follow the formalism outlined by Hughes (2002). 



1 University of Melbourne, Parkville, Victoria, Australia 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St.. 

3 



Cambridge, MA 02138 



http: / / www.ligo.caltech.edu/advLIGO / 



2 







i i I i i 


o 




■ '"■ '-^&4 Zq= 


M se e d = 


260M so|Qr : 


CM 


. \ ^ q =5 M seed = 20M solar : 


... . v w 






o 

CM 


:. '. '•>'.'. A 

r \> - - 






, . i , ; 


o 

to 

CM 

o 

CM 
O 


E i ■ ■ ■■■ ■ .i . ; 


10 




100 




10 100 


^rd,obs 


(Hz) 






f rd,obs ( Hz ) 



Fig. 1. — Comparison between the ring-down strain signals and the LIGO-II sensitivity curve as a function of observed frequency, 
/rd obs- The points correspond to hierarchical coalescence of 10 MBH seeds to a single MBH at a redshift Zq. Two cases are shown, having 
MBH seeds of 20Mq and 260Mq, coalescing by z q = 5 and z q = 1, respectively. Each coalescence event involves a different randomly-chosen 
combination of the parameters 9, <j>, ip, A and L ■ h. The calculations assume e r( j = 0.03 and a = 1.0. 



The dimcnsionlcss strains for the two polarizations of the 
gravitational waveform during the ring down may be writ- 
ten as 

h+,-(t) = A+- CXp (-Ttfrd.obst/Q) cos (27r/ rdi0bs i) , (1) 

where Q = it / r d,obsT bs is the quality factor, r b s is the 
observed decay time and t is the observed time following 
the merger. The observed frequency, / r d,obs, equals the 
intrinsic frequency, / r d, divided by (1 + z). Assuming that 
the ratio of the strain amplitudes for the two polarizations 
follows that of the inspiral phase, 



A+ = Aing 1 + (L- n) 2 



; A, 



-2Arin S {L-n), (2) 



where L ■ n is the dot product of the unit vectors pointing 
along the binary orbital angular momentum and the line 
of sight to the source. The overall strain amplitude _4 rc i is 
found by requiring that the binary radiates a fraction e rd 
of its energy (e.g. Fryer, Holz & Hughes 2002). If a black 
hole of mass AMbh merges with a black hole of mass Mbh 
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(3) 

Note that in this expression D — D^/(l + z) is the co- 
moving distance, where £>l is the luminosity distance 
(e.g. Hogg 1999). The energy flux is proportional to 

(/rd,obs) 2 ^ ing « [(/rd) 2 (l + Z)- 2 ] D- 2 CX Dl 2 . ^ 

Distortions of the Kerr metric during the ringing phase 
may be decomposed into spheroidal modes with spher- 
ical harmonic like indicies / and m (Fryer, Holtz & 
Hughes 2002). The quadrupole (I — 2) moments are ex- 
pected to dominate ringdown spectrum. Binary coales- 
cence excites a bar-like deformation of the event horizon 
(Hughes 2002). As a result, the m = 2 mode dominates 
over the m — mode when the ringdown forms the final 
stage of a binary merger. Fits to numerical simulation 
show that for l,m = 2 (Leaver 1985; Echeverria 1989; 
Fryer, Holtz & Hughes 2002) 
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(4) 
(5) 



where a is the dimensionless black-hole spin, taken to have 
representative values of 1, 0.5 or in the numerical results 
presented later. Hughes & Blandford (2002) have found 
that equal mass mergers may have rapidly rotating rem- 
nants (a ~ 1), while mass ratios of 0.5 lead to typical 
spins of a ~ 0.5. We note that in addition to setting the 
oscillation frequency, the spin governs the value of Q and 
hence the number of oscillations. Ringdown waveforms 
with more oscillations will be easier to detect. 
The observed strain may be written as 

H (t) = F+(0, 0, 1>)h+(t) + F x (6, <f>, VO/ix (t), (6) 

where F + (9 7 (f) 7 ip) and F + (6,(j),ip) are the detector re- 
sponse functions for a source with a sky position (6, </>) and 
polarization axes rotated at an angle ip (e.g. Thorne 1987). 
Since the ring-down waves lie in a narrow frequency 
range, the signal-to-noise ratio p for a detection uti- 
lizing matched-filtering techniques may be evaluated as 
(Hughes 2002) 

2 _2f °°dtH(t) 2 
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(7) 



where 5h(/ r d,obs) is the spectral density of detector noise 
(in Hz" 1 ). 

3. MBH MERGERS AND THE EVOLUTION OF THE SMBH 
POPULATION 



In the redshift interval between z q and z„ 



dzq, the 



accretion due to quasars implies that the co-moving mass 
density in SMBHs is 

f°° Lu i dt 
dp hh (z q )=dz q dL B — ^-^(LB,z q ) — , (8) 
Jo ec ®z q 

where Lboi is the bolometric luminosity of a quasar with 
a B-band luminosity Lb, e is the conversion efficiency 
of accreted mass to electromagnetic radiation (we adopt 
e = 0.1; Yu & Trcmainc 2001), and *(L B ,z q ) is the co- 
moving density of quasars per unit B-band luminosity at 
z q , for which we adopt the parametric form derived from 
the 2dF survey (Boyle et al. 2000). By integrating equa- 
tion (8) over redshift, one may obtain the SMBH mass- 
density accreted by a redshift z q . The local (z q — 0) value 
for the co-moving mass density in SMBHs accreted dur- 
ing the quasar activity is p bb w 2.7 x 10 5 M Q Mpc~ 3 . This 
number is similar to estimates of the local SMBH mass 
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density obtained from the inventory of quiescent galac- 
tic SMBHs (e.g. Yu & Tremaine 2002; Aller & Rich- 
stone 2002). 

If some fraction of the SMBH mass is assembled through 
coalescence during the hierarchical build-up of galaxies, 
then it is natural to expect MBH remnants to contribute 
to the SMBH mass budget as they are more massive than 
normal stellar remnants and should sink to the centers of 
galaxies by dynamical friction (Madau & Rees 2001). The 
initial coalescence of the MBHs would be detectable by Ad- 
vanced LIGO, and the observed event rate (or lack thereof) 
could shed light on the possibility that MBH coalescence 
events produced the early seeds of SMBHs. Suppose that 
a fraction /mbh of the mass in SMBHs at any redshift was 
added through coalescence of MBHs, with the remainder 
added through gas accretion. Since the comoving mass 
density of accreted gas mass for quasars is comparable to 
the final SMBH mass density in the local universe, we ex- 
pect /mbh "C 1 independent of other considerations. 

While the accreted mass is being added at redshift z q , 
the coalescence of the MBHs incorporated into the SMBH 
at z q would have stretched over most of the age of the 
universe t ag e(z q ), corresponding to the assembly time of 
the SMBH. The MBHs are brought together through hier- 
archical merging of the mini-halos that form the building 
blocks for the massive dark matter halo hosting any par- 
ticular SMBH. There are several issues to consider, includ- 
ing the number and clustering of MBHs inside mini halos, 
the merger rate of mini-halos into proto-galaxies and more 
massive hosts, and the subsequent dynamical behavior of 
the MBHs within the resulting gaseous and stellar environ- 
ments. Rather than model the coalescence history using a 
dark-matter halo merger tree and the poorly understood 
dynamics of the MBH population, we instead assume that 
the MBHs undergo random hierarchical coalescences with 
the MBH mergers distributed linearly in time. At any 
point in the coalescence history, two MBHs that coalesce 
are chosen at random from the MBH population. The next 
merger is then drawn from the revised population which 
contains one fewer MBH in total, and one additional MBH 
of a larger mass, and so on. There are a total of N — 1 
coalescence events leading from the N initial seed MBHs 
to the single final MBH of mass AMmbh- One third of 
these coalescence events are between MBHs of the initial 
seed mass. 

Figure 1 shows the signal amplitude H<\ = 
\J 2 J °° dtH(t) 2 for the hierarchical coalescence of 10 3 

MBH seeds to a single MBH at a redshift z q . Two cases 
are illustrated, having MBH seeds of 20M Q and 26OM 



coalescing by z q 



5 and z Q 



1, respectively. Each 



individual coalescence event was combined with a differ- 
ent randomly chosen combination of the parameters 9, 
<j), ip, A and L ■ n. The calculations assume e r d ~ 0.03 
(Hughes 2002), and a = 1.0. The values of Hd in Figure 1 
lie along two branches. The lower branch originates from 
the early stage of the merger tree, where the seed MBHs 
provide most of the coalescence events. As the typical 
mass of the coalescing MBHs grows at lower redshifts, 
the signal amplitude increases while its frequency de- 
creases. The figure also shows ^/5h(/ r d,obs), the planned 
sensitivity curve for the two LIGO-II 4km interferome- 



ters (Gustafson, Shoemaker, Strain & Weiss 1999). The 
signal-to-noise ratio p is the ratio between the points and 
the sensitivity curve [see equation (7)]. From an ensemble 
of such detection histories, we may therefore compute the 
fraction of mergers /detect (-Zqj Zcoai) that are detectable by 
LIGO-II as a function of the quasar redshift, z q , and the 
coalescence redshift, z coa i. 

The rate at which coalesced MBH mass is added to the 
SMBH population is given by 

rfpMBH /mbh dpbYi 



dz a 



1 



, ■ (9) 

^q J- — /MBH «Z q 

with the remainder added via gas accretion during the 
active quasar phase. The number density increment of 
consumed MBHs can be derived from the mass density 
increment, dn M Bn(z q ) = dpMBn(z q )/M M BH- 

The redshift distribution of the rate of detectable coa- 
lescences (in the observer frame) is 
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(10) 



-^coal , 

where the factor 1/(1 + z coa i) is due to cosmological time 
dilation, the derivative dt/dz coa \ enforces a linear distri- 
bution of coalescence events in time, (dV/dQdZq) is the 
co- moving volume per unit solid angle per unit redshift, 
and 0(a;) is the Heaviside step function. 

Figure 2 shows d 2 Nd ct /dtd\n(l + z coa i) in units of 
10 /mbh yr _1 - The plotted curves correspond to a signal- 
to-noise ratio of 5, assuming a MBH initial mass of 2OM 
(left panels) and 26OM (right panels). The calculations 
were performed for values of a — 1.0 and e r d = 0.01, 
0.03 and 0.10. A fraction as small as /mbh ~ 10~ 5 
yields ~ 0.1-10 detectable coalescence events at redshifts 
2 < z q < 5 (where the quasar population is substantial 
and the observed frequency does not redshift out of the 
Advanced LIGO band above 10 Hz). For comparison, 
the light lines show the redshift distribution of the en- 
tire coalescence rate (/detect = 1) in the observer frame, 
d 2 N mlg /dtdln(l + z coal ). We find that Advanced LIGO 
will be able to detect most of the low redshift coalescences, 
while the detectable fraction declines at early cosmic times 
due to the redshift in frequency. 

Values of the total event rate {dN^et/dt) are presented 
in Table 1 in units of 10 5 /mbh yr _1 . In addition to the 
examples plotted in Figure 2, the table lists results for 
BH spins of a = 0.0 and a = 0.5. The rates are lower for 
smaller MBH spins and for larger MBH masses, due to the 
lower emission frequency and smaller numbers of MBHs. 

4. DISCUSSION 

We have estimated the number of seed MBHs that must 
have coalesced at high redshifts in order to contribute a 
fraction /mbh of the mass density in quasar SMBHs. We 
then found the corresponding detection rate of ring-down 
gravitational waves for the Advanced LIGO interferome- 
ters. If the initial mass of MBH remnants was 2OM , then 
1 event per year would imply a coalesced mass fraction of 
/mbh ~ 1CT 6 . On the other hand, if the MBH seed mass 
was 26OM , then 1 event per year would imply a coalesced 
mass fraction of /mbh ~ 10 -4 . Advanced LIGO will be 
primarily sensitive to MBH coalescence at z ~ 2 — 3, while 
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Fig. 2. — The rate of detectable mergers per ln(l + z) in units of 10 5 /mbh yr 1 ■ Curves are shown for a signal-to-noise ratio of 5, and 
initial MBH mass of 2OM (left panels) and 260Mq (right panels). Each panel shows three curves for e rd = 0.01 (solid), 0.03 (dashed) and 
0.10 (dot-dashed). The calculations assumed a = 1.0. For comparison, we also show the total number of mergers per year in the observer 
frame (light lines). 



Table 1 

Detection rate of coalescence events in units of 10 b /mbhYR _1 . The third column shows the merger rate 

assuming that all events can be detected. 

[l(r b yr// M BH] x dN dct /dt 



seed 


a 


[10- b //mbh] x dN^g/dt 


e rd = 0.01 


e rd = 0.03 


e rd = 0.10 


2OM 


1.0 


35.6 


1.9 


5.7 


13.3 


20M Q 


0.5 


35.6 


1.0 


3.4 


9.8 


2OM 


0.0 


35.6 


0.67 


2.3 


7.5 


26OM 


1.0 


2.7 


0.02 


0.07 


0.15 


26OM 


0.5 


2.7 


0.003 


0.01 


0.03 


26OM 


0.0 


2.7 


0.001 


0.003 


0.007 



our calculation implies that most of the MBHs for which 
coalescence is possible may have coalesced earlier. How- 
ever, values of /mbh ~ 10~ 4 — 1CP 6 correspond to only a 
few events per typical galactic nucleus per Hubble time, a 
very small coalescence rate that is comparable to the rate 
by which massive galaxies merge. 

It is possible that not all MBH coalescences contribute 
to the assembly of the SMBH mass. In particular, MBHs 
that are born into tight binaries might coalesce, even if 
they do not migrate into dense stellar systems or MBH 
clusters. If MBH binaries did coalesce without adding 



mass to the assembled SMBHs, then the observed event 
rate would only provide an upper limit on /mbh- MBHs 
that coalesced directly with SMBHs without first under- 
going coalescence with other MBHs would be invisible to 
Advanced LIGO, but might be detected by the planned 
space mission LISA 4 (Wyithe & Loeb 2003). 
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AST-0204514 for AL. 
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